Orbital-selective Mott Transitions in a Doped Two-band Hubbard Model 
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We extend previous studies on orbital-selective Mott transitions in the paramagnetic state of the 
half-filled degenerate two-band Hubbard model to the general doped case, using a high-precision 
quantum Monte Carlo dynamical mean-field theory solver. For sufficiently strong interactions, 
orbital-selective Mott transitions as a function of total band filling are clearly visible in the band- 
specific fillings, quasiparticle weights, double occupancies, and spectra. The results are contrasted 
with those of single-band models for similar correlation strengths. 
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I. INTRODUCTION 

Starting with the early theoretical work by Hubbard^ 
and Brinkman and Rice,— and experimental work on tran- 
sition metal oxides,^ the Mott-Hubbard metal-insulator 
transition (MIT) has been a subject of great interest in 
condensed matter physics for decades. Much of its es- 
sential physics is captured already by the single-band 
Hubbard model^'^'^ as has been shown in numerous 
studies using the dynamical mean-field theory (DMFT) 
frameworkii^ Multi-band extensions of the Hubbard 
model allow for a more realistic description of the MIT 
and of other strong-correlation phenomena and have, 
in general, much richer phase diagrams. Multi-orbital 
Mott transitions have been explored (within DMFT) in 
a doped two-band model with equivalent orbitals already 
more than 10 years a.go^ In this case, the band degener- 
acy increases the number of Mott lobes from 1 to 3, but 
has no fundamental impact on the Mott physics near half 
filling (in contrast, e.g., to ferromagnctism for which the 
inter-orbital Hund exchange is essentialiSiiiii^iHiii). 

More recently, it has been realized that orbital de- 
generacy can change the character of Mott transitions 
in a fundamental way: namely, in the case of multi- 
ple inequivalent orbitals. In the electronic context, such 
inequivalence can arise naturally by orbital-dependent 
hopping amplitudes, associated, e.g., with in-plane ver- 
sus out-of-plane t2g orbitals in layered ruthenatesj ^^i^^ 
It has been suggested by Anisimov et alJ^ that such 
systems should undergo a sequence of orbital-selective 
Mott transitions (OSMTs) with increasing interaction 
strength, which would then explain the peculiar phase 
diagram of Ca,2-xSTxRu04r^^^ Clearly, such a scenario 
- with coexisting itinerant and localized valence electrons 
in the intermediate orbital-selective Mott phase - is fun- 
damentally different from a conventional simultaneous 
Mott transition of all valence electrons. 

Within the last few years, it has been estab- 
lished that OSMTs occur in multi-band Hubbard 
models with inequivalent orbitals under quite general 
circumstancesJ^ii^i^i^i22ii^— A (nearly) mini- 
mal Hamiltonian for orbital-selective Mott behavior is 
the two-band Hubbard model with band-specific hopping 
amplitudes tm (for orbital index m G {1,2}) and Ising 



type Hund rule couplings (parametrized by Jz), 
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Here, a G {t, 1} denotes the spin; i and j label lattice 
sites; fiima = c\^„c^m(j- contrast to the hopping, the 
intra-orbital Hubbard interaction U is assumed orbital 
independent. The intra- and interorbital interactions are 
related by [/ — U' + 2Jz. In the following, we will refer 
to Eq. (HD as the J^-model,^i^i^ 

A general discussion of the Hund exchange would have 
to include spin-flip and pair hopping terms. 



However, such terms arc not essential for OSMTs (Ref. 
[28l ) and change their character only in the SU(2) sym- 
metric limit Jz = J± = J (i.e., in the Heisenberg limit). 
In particular, it was shown for the half-filled casa^ that 
the itinerant electron species in the orbital-selective Mott 
phase is a non-Fermi-liquid for all J± < Jz and a Fermi 
liquid only at J± = Jz- In this sense, Jz model ([T]) rep- 
resents the generic class of anisotropic Hund couplings. 

The nature of OSMTs and their essential requirements 
are well understood for particlc-holc symmetric situa- 
tions. This symmetry is destroyed upon doping (i.e., by 
variation of the chemical potential fi) and/or by adding 
an orbital dependent field, i.e., by crystal field splitting. 
Koga et al. studied effects of weak hole doping on OSMTs 
in a two-band Hubbard modeli^ Crystal field effects have 
been investigated in two-orbital2S and three-orbital^ sys- 
tems. In a ferromagnetic slave-boson mean-field theory 
calculation, Riiegg et ali^ obtained a phase diagram for 
the two-band Hubbard model with crystal field splitting, 
which includes ferromagnetic as well as band insulating 
regimes. 

In this paper, we extend high-precision^ quantum 
Monte Carlo (QMC) calculations22i30i36 ^jthin DMFT 
for two-band model ([1]) to general doping levels. Of cen- 
tral interest is the phase diagram, which is derived from 
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orbital-dependent fillings, double occupancies, and quasi- 
particle weights. At the same time, these observables as 
well as spectra are used for characterizing the orbital- 
selective physics, i.e., the effects of orbital inequivalence. 
Other results of interest in this paper include unexpected 
behavior in the intermediate interaction regime, where 
the double occupancy for the wide band becomes nearly 
independent of the interaction. We also show that the 
entire doping range < n < 4 is subdivided into three 
different regimes by the transport behavior of the narrow 
band, and that the two-band spectra differ from analo- 
gous single-band results in that the additional scattering 
channels introduced by the interorbital couplings tend 
to smear out the spectrum. Finally, we will discuss the 
relation between orbital-selective Mott physics and non- 
Fermi- liquid behavior, i.e., whether these phenomena al- 
ways occur simultaneously (in the model). Here, as 
in previous work on this subject, we concentrate on the 
effect of correlations on the thermodynamic and spectral 
properties of the system in the paramagnetic phase, i.e., 
excluding magnetism. 

The structure of this paper is as follows: In Section IlSI 
we briefly discuss the methods and model parameter val- 
ues relevant to our work. Section II Bl gives information 
on the observables studied in our paper. Static and dy- 
namic observables are treated separately (in Sees. |ll] and 
mil respectively). The static observables discussed in Sec. 
im include band-resolved particle numbers fSec. Ill A|) . in- 
traorbital double-occupancies (Sec. IIIBj) . and the phase 
diagram as a function of interaction and total filling (Sec. 
Ill C[) . As dynamic observables, we study band-resolved 
spectral functions (Sec. IlllXj) . the spectral weight at the 
Fermi edge (Scc. lIIIBl) . and the Matsubara self-energy as 
well as quasiparticle weights (Sec. IIII Cl) . Finally, in Sec. 
IIVI we summarize the results and formulate our conclu- 
sion. 



A. Methods 

In the following, we consider J^-model ([T|) for two 
bands with a bandwidth ratio of 2. For the narrow band, 
we assume a semi-elliptic "Bethe" non-interacting DOS 
with a fixed hopping amplitude tn = 0.5 and, hence, a 
full bandwidth Wn = 2. Similarly, we assume a semi- 
elliptic DOS with iw = 1 and W^, = 4 for the wide band. 
The intraband on-site interaction U is chosen as the pri- 
mary (variable) interaction parameter; it acts equally in 
both orbitals. The interorbital interactions are scaled 
as Jz = U/A and U' = U/2, in line with most earlier 
studies jiSii^iSSii^i This work is restricted to the param- 
agnetic case; spin indices are omitted when appropriate. 
We choose a fixed temperature T = 1/40, which is of the 
order of the critical temperature in the undoped casei^ 

Within DMFT, two-band model ^ is mapped to a 
two-orbital single-impurity Anderson model which has 
to be solved self-consistently4 In this work, the im- 
purity model is solved using the Hirsch-Fyc quantum 



Monte Carlo (HF-QMC) algorithm,^i2S This method 
discretizes the imaginary-time path integral expression 
for the Green function into A time slices of uniform width 
At = /3/A, where /? = 1/T (for fee = 1); a Hubbard- 
Stratonovich (HS) transformation replaces the electron- 
electron interaction at each time step by a binary auxil- 
iary field which is sampled using standard Markov Monte 
Carlo techniques. The results are exact in the combined 
limit of vanishing discretization At — » and a large 
number of Monte Carlo update sweeps. In this work, 
the systematic errors associated with finite discretization 
(At = 0.4, unless noted otherwise) are minimized by 
supplementing the discrete QMC Green function with a 
high-frequency expansion of the self-energy i'^^i^" The high 
precision of this method has proved essential in detecting 
both OSMTs in the undoped cascj^ if required, extrap- 
olations At can be used in order to achieve extreme 
precision^ and efficiencyi^ In order to compute spectral 
functions, the imaginary-time Green functions resulting 
from QMC are analytically continued to the real axis us- 
ing a standard maximum-entropy method (MEM) 4^ 

We have verified that the remaining discretization bias 
(at At — 0.4) is generally small (much smaller than in 
preceding QMC implementations^) , up to slight shifts in 
the critical interactions, by performing additional simu- 
lations with 0.25 < At < 0.5. In particular, we did not 
find significant effects on the single-particle spectra: their 
general shapes as well as peak positions and weights were 
unchanged. 



B. Observables 

Using the QMC-DMFT methods, described above, 
we investigate several observables in this paper: band- 
specific particle numbers and intraorbital double occu- 
pancies, the density of states, in particular its value at 
the Fermi level, and the quasiparticle weight. We now 
briefly introduce these quantities, which are of interest by 
themselves and are used in order to construct the phase 
diagram. 

To start with the general properties of J^-Hamiltonian 
([T]): this model is formulated in a particle-hole symmetric 
way (in particular, no crystal field splitting is included), 
so that results for a total filling of n > 2 can be easily 
mapped to results for n < 2. As a result of this well- 
known symmetry, particle- and hole-dopings are equiva- 
lent and, accordingly, the word "doping" in the following 
refers to both. Moreover, the phase diagrams, to be cal- 
culated below, will be mirror symmetric with respect to 
the half-filled case, where n = 2. 

The study of orbital-specific particle numbers and in- 
traorbital double occupancies is particularly interesting 
for parameter values, for which the orbitals behave phys- 
ically differently, i.e., in the orbital-selective Mott phase 
(OSMP). In this case, the narrow band is in the insulat- 
ing state and displays a gap, while the wide band is in a 
non-Fermi-liquid metallic state. The insulating nature of 
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the narrow band physically suggests that both the densi- 
ties and the double occupancies of this band may initially 
remain pinned upon doping, and in fact one of the goals 
of our investigation of these observables is to determine, 
whether and to what extent (depending upon the inter- 
action strength U) this is true. 

Clearly, in addition to the transport properties (e.g., 
the electrical conductivity), which are not explicitly in- 
vestigated in this paper, one of the most fundamental ob- 
servables for detecting metal-insulator transitions is the 
density of states (DOS), i.e., the local spectral function. 
It is particularly interesting to contrast the results for 
the DOS in the two-band calculations with single-band 
results for comparable parameter sets, whenever appro- 
priate. 

Another important indicator for the onset of insulating 
behavior is the quasiparticle weight Z, which is mathe- 
matically determined by the self-energy S(i^). Upon ap- 
proaching the insulating phase of a single-band model 
from the metallic side, the quasiparticle weight vanishes 
(at least in the ground state) at the metal-insulator tran- 
sition. On account of the relation Z = m/m*, the van- 
ishing of the quasiparticle weight is equivalent to the di- 
vergence of the effective mass m*. It is important to 
note that the physical interpretation of Z and m* as 
the "quasiparticle weight" and the "effective mass", re- 
spectively, is meaningful only in a Fermi-liquid regime. 
Within a Fermi-liquid phase, the quasiparticle weight is 
well approximated in an imaginary-time calculation by 
its imaginary-time analog: 



m 

771* 



^ 1- 



duj 
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ttT 
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Here uJi ~ ttT is the first Matsubara frequency. Both ex- 
pressions for the quasiparticle weight clearly recover the 
non-interacting limit Z — > 1 at arbitrary temperatures 
for S — > 0. In addition, the imaginary-time secant ap- 
proximation is exact in a Fermi-liquid phase in the limit 
T — > 0. However, both expressions are demonstrably dis- 
tinct, e.g., in the insulating phase at T = 0, since the 
real-frequency definition of Z vanishes, whereas estimate 
([2]) remains finite. This illustrates that estimate ([2]) will 
in general be inaccurate or even invalid in non-Fermi- 
liquid regimes, such as in particular the OSMP, where 
the narrow band is insulating, while the wide band is 
a non- Fermi-liquid type metal (see). In any case, the 
discrete imaginary-time expression in Eq. Q is an es- 
tablished measure of the low-frequency behavior of the 
self-energy. 



II. STATIC OBSERVABLES AS A FUNCTION 
OF DOPING 

The total electron concentration in a two-band model 
is composed of contributions from the two orbitals: n = 
"w + rin, where riw = I]o.(^twcr) and Z]cr(^mo-) 
represent the particle filling of the wide and the narrow 
band, respectively. 

In two-band system ([l} with total half filling, due to 
particle-hole symmetry, either band is exactly half-filled: 
n-w = n-n = 0.5. Let us briefly review the central results 
obtained earlier for this special casei^^ with increasing 
interaction, and starting from the non-interacting limit 
(and for low T — 1/40), both bands undergo consecutive 
Mott transitions. For weak interaction {U < Ud ~ 2.1), 
both bands remain itinerant. For strong interaction 
{U > Uc2 ~ 2.6), both bands become insulating, since 
the energy difference between the atomic levels becomes 
comparable to the interaction U. In the intermediate re- 
gion, the wide band is metallic and the narrow band is 
insulating. Slight variations in temperature do not signif- 
icantly change the results quantitativelyi^ This picture 
agrees with arguments in Ref. [13, relating the interacting 
two-band model to two distinct bands, which separately 
undergo the Mott transition according to their individual 
bandwidth. 

In the following, we first discuss our results for the 
band-specific particle numbers and double occupancies, 
which are then compared to predictions from a simple 
rigid-band model and from perturbation theory. In com- 
paring the QMC results to simple model calculations, 
the doping dependence turns out to be of great inter- 
est. We also present a phase diagram for the doped 
two-band model, with particular emphasis on transport 
characteristics. Accordingly, we distinguish insulating, 
orbital-selective Mott and metallic phases. 

A. Band-specific particle numbers 

We start the presentation of our results with the band- 
specific fillings, which are plotted in Fig. [1] for different 
interaction strengths 1.8 < U < 2.8. For small deviations 
from half filling and U > 2, only the wide band accounts 
for the doping, while the narrow band remains (to a very 
good approximation) half-filled. As a consequence, the 
relation ~ n — 1 holds approximately in this low- 
doping regime. As one expects, the pinning of the parti- 
cle filling of the narrow band at its half-filled value in the 
low-doping regime becomes more extended for increasing 
U. These plateaus are bounded by well-localized kinks, 
beyond which the curves ri,i(n) show rapid, nearly lin- 
ear, increase: additional electrons are predominantly al- 
located to the narrow band. Hence, the slope of n^{n) 
drops significantly. The corresponding phase boundary 
will be discussed below, in Sec. Ill CI Note that the curves 
in Fig. [1] for the narrow and the wide bands cross in the 
OSMP regime near n ^ 2.4, with details depending upon 
U. 
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FIG. 1: (Color online) Band-specific fillings for the two-band 
model versus total filling for a range of interactions U. Open 
and solid symbols represent the narrow and the wide band, 
respectively. 



B. Double occupancies 

In Fig. [2^, we show QMC results for the intraorbital 
double occupancies (nimt^-imi) (™ G {w,n}) of the wide 
band (upper curves) and the narrow band (lower curves). 
In the weakly doped regime {n < 2.1) of the OSMP 
(with U ^ 2.2), the intraorbital double occupancy for 
the narrow band - which, as was demonstrated above, 
is effectively singly occupied - is strongly suppressed. 
In this regime {n < 2.1), the double occupancy of the 
wide band rapidly decreases with increasing interaction 
due to Coulomb repulsion. In the stronger doped regime 
(n ^ 2.1) the situation reverses: The double occupancy 
is supressed in the narrow band with increasing interac- 
tion, whereas the intraorbital double occupancy for the 
wide band becomes nearly U independent. As we will see 
below, both bands are metallic in this regime (n ^ 2.1) 
and show a significant DOS in the vicinity of the Fermi 
edge, the DOS of the narrow band obviously being larger 
than that of the wide band. As a consequence, both the 
average filling and the concentration of doubly occupied 
orbitals increases more rapidly in the narrow band than 
in the wide band for n > 2.1. Hence, as for the particle 
numbers, a crossing of the curves for the wide and nar- 
row bands occurs at n ~ 2.5, so that the concentration 
of intraorbital double occupancies is larger for the nar- 
row than for the wide band for n > 2.5. Moreover, these 
curves show a very weak [/-dependence in this regime, 
which suggests that very similar behavior is to be ex- 
pected from weak-coupling perturbation theory. We will 
come back to this issue below. The crossing and the 
virtual U independence can be clearly seen in Fig. [^b, 
which shows the QMC results for the double occupan- 
cies in the entire regime 2 < n < 4; for comparison, the 
corresponding weak-coupling (Hartree-Fock) and nonin- 
teracting results are also shown. 



a) 0.4 




n (total filling) 

FIG. 2: (Color online) Intraorbital double occupancy for var- 
ious [/-values as a function of doping for the two-band model. 
Panel (a) shows for dopings near half filling and panel (6) in 
the entire doping range 2 < n < 4 that the intraorbital double 
occupancy for the wide band becomes nearly independent of 
the interaction for intermediate U and all fillings n > 2.2. In 
(6), results for the noninteracting case ([/ = 0, black line) and 
within Hartree-Fock approximation (for U = 2.4, long-dashed 
line) are included for comparison. 



The primary mechanism governing the band-specific 
fillings and the intraorbital double occupancies in the 
low-doping regime can be understood on the basis of a 
rigid-band model with two bands and intraorbital inter- 
action only. A comparison between the band-specific fill- 
ings within this rigid-band model [sketched in Fig. [3] with 
the QMC results reveals good agreement for n < 2.1, so 
that the pinning of the narrow-band filling can be ex- 
plained as resulting from the single-particle gap. If one 
further assumes that the narrow-band double occupancy 
is a function of the narrow-band filling and of the inter- 
actions, it has also to remain pinned in the same regime 
- in full agreement with the QMC results. The validity 
of the rigid-band model rapidly breaks down as the total 
filling is increased beyond n > 2.1, where the actual par- 
ticle numbers of the narrow (wide) band increase more 
rapidly (slowly) than in the rigid-band model. The de- 
viations between the simple rigid-band assumption and 
the true results in the metallic phase (n > 2.1) are due 
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FIG. 3: (Color online) (a) DOS of a rigid-band model with a 
semi-elliptic DOS and bandwidths Wn = 2 and Ww = 4, rep- 
resentative of interactions U ~ 2.4. (b) Band-specific fillings 
within this model, compared to QMC results for U = 2.4 and 
T = 1/40. 



to correlations, in particular to the formation of Kondo 
resonances in both bands. 

We now discuss and interpret the remarkably weak U- 
dependcncc of the intraorbital double occupancies of the 
narrow and the wide band for n ^ 2.5. Before compar- 
ing to perturbative results, we illustrate the weak depen- 
dence of the double occupancies on U in Fig. [H which 
covers the interaction regime < U < 7 for a total den- 
sity n = 3. Figure |4] indeed shows a near-constancy of 
the double occupancies for such interaction values. On 
the basis of the virtual {/-independence, in particular for 
2 ^ t/ ^ 4, one expects good agreement between QMC- 
results and perturbation theory in this regime. As can 
be seen from Fig. [2)d, where corresponding Hartree-Fock 
results are also plotted, this is indeed what one finds. 
Figure [21d shows that the wide band agrees with per- 
turbation theory better than the narrow band, at least 
for n < 3, as one expects, since the effective interaction 
strength (compared to the band width) is larger for the 
narrow than for the wide band. This is also confirmed 
by Fig. [5l which shows the quantum fluctuations around 
the Hartree-Fock solution, which (for n < 2.7) are indeed 
largest for the narrow band. The near-perfect agreement 
between the narrow-band QMC results and perturbation 
theory for n > 3.2 follows immediately from the near ab- 
sence of singly occupied or empty orbitals in the narrow 
band in that regime. For completeness we should add 
that the weak U -dependence of the concentrations of in- 
traorbital double occupancies clearly cannot persist for 
arbitrarily large U values. From strong-coupling pertur- 
bation theory one expects that, for sufficiently large U, 
the double occupancies of both bands converge towards a 
common strong-coupling value 0.5. 

In between such a strong-coupling phase and the 
metallic weak-coupling phase, one would expect a metal- 



FIG. 4: (Color online) Intraorbital double occupancy ver- 
sus U for density n = 3 (squares). In the interaction range 
2.5 ^ U 5, the intraorbital double occupancy is only very 
weakly dependent on U. Density-dependent upper and lower 
bounds on the double occupancies are obtained in the limits 
of weak (dashed/solid lines) and strong (dotted lines) corre- 
lation, respectively. Inset: quasiparticle weight Z (cf. Sec. 

img. 



insulator transition. As mentioned in the introduction, 
such additional Mott lobes (at quarter and three quarter 
filling) had been found in an early HF-QMC study of a 
two-band Hubbard model, even with critical interactions 
similar to those at half fiUingi^ However, this study con- 
sidered the case of SU(4) symmetric interactions {U = U' 
and J = 0) and equal bandwidth. The first difference 
is particularly relevant, since the Heisenberg exchange 
strongly suppresses Mott phases at (three) quarter filling 
in terms of the critical interactions: the phase boundaries 
are shifted to much larger Uc- This is easily seen by com- 
paring the cost of a charge fluctuation in the atomic limit 
{jU at n = 1 or n = 3, compared to |J7 at ti = 2; both 
for J = U/A,U' = U/2) and has also been found numer- 
ically for SU(2) symmetric interactions4^ In the present 
case, we expect the system to remain metallic (at n = 3) 
up to Uc ^ 20, i.e., beyond the interaction range easily 
accessible using QMC calculations. In addition, the crit- 
ical temperatures may be considerably lower; however, 
sharp crossover lines should persist at least up to the 
temperatures studied in this paper 4^1^ 

From the data of Fig. 21 we can, indeed, exclude a 
phase transition at n = 3 for U 7: both the intraor- 
bital double occupancies (squares, main panel) and the 
quasiparticle weights Z (inset) are completely smooth as 
a function of the interaction [/, i.e., do not show any 
kinks or hysteresis effects that would be characteristic of 
an MIT. Note that both the double occupancy and Z 
are larger for the narrow band due to its larger filling 
fraction. 
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FIG. 5: (Color online) The t|-density-density correlation 
function for total fillings n > 2.2, corresponding to quan- 
tum fluctuations in the intraorbital double occupancy around 
the uncorrelated limit. The correlations are generally nega- 
tive, since an enhanced |-density on a site implies a decreased 
j-density. At fixed interaction U, the amplitude of the corre- 
lation function decreases monotonically with filling, since the 
available phase space for density fluctuations decreases with 
increasing n > 2. 

C. Phase diagram 

Next we discuss the phase diagram as a function of the 
interaction strength U and the total fiUing n with partic- 
ular emphasis on transport-related properties. To do so, 
we distinguish an insulating state of the two-band sys- 
tem from an orbital- selective Mott phase, in which one 
band is metallic and the other insulating, and a purely 
metallic phase, in which neither of the bands is insulat- 
ing. As is well-known, all three phases occur already at 
half filling, where the system is in the metallic, orbital- 
selective Mott or insulating state for U < Ud — 2.1, 
Uci <U< Uc2 — 2.6, or Uc2 < U, respectively. 

For the doped system {n ^ 2) it is already clear from 
our investigation of particle numbers and of double oc- 
cupancies for the narrow band at low doping, which ap- 
peared to be pinned at the half-filled state for moderate- 
to-large [/-values, that this band is insulating in this 
regime. Since the particle filling of the wide band changes 
substantially with total density (according to — Ji — 1) 
for these {U, rt)-values, one expects the wide band to be 
metallic in this part of the phase diagram. Accordingly, 
one can take the pinning of the narrow-band particle 
numbers as a criterion for the occurrence of an OSMP. 
Alternatively, at larger doping, where the fillings of both 
bands change upon changing the total density, it is clear 
that both bands have finite compressibility and, hence, 
are metallic. Finally, at non-integer total filling in the 
paramagnetic phase, one does not expect the occurrence 
of a rigorously insulating state for both bands simulta- 
neously, since at least one of the bands must have fi- 
nite compressibility in that region. Thus, to summarize. 




1.8 2 2.2 2.4 2.6 2.8 3 

U (interaction) 



FIG. 6: (Color online) Orbital-selective Mott phase diagram 
of the two-band model for particle densities 1.85 < n < 2.15. 
For interaction U < Ud, both bands are metallic at all den- 
sities. The insulating phase at half filling (thick black line) 
only exists for U > Uc2- 



we found a simple criterion for the determination of the 
paramagnetic phase diagram as a function of interaction 
and doping: If the system shows pinning for the narrow 
band away from half filling, it is in an OSMP, if not, then 
it is metallic. 

It is clear that this criterion for the determination of 
the phase diagram, which is essentially based on the be- 
havior of the compressibility as a function of U and n, 
should be cross-checked with alternative criteria, based 
on dynamical quantities, such as. e.g., the finitcness of 
the DOS at the Fermi level for the wide and the nar- 
row bands. Results for such dynamical quantities, which 
are presented below, turn out to be fully consistent with 
the present "static" criteria and lead to the same phase 
diagram. 

The phase diagram, which is based on the above- 
mentioned criteria concerning the band-specific filling 
factors nn(n) and ny^{n) as a function of the total density, 
is presented in Fig. [5] for the moderately strong interac- 
tions U < 3, near the OSMP previously established at 
half filling. One finds that the purely insulating state, 
which occurs only at half filling as explained above, is 
embedded in an OSMP, which reaches from U ~ 2.1 to 
U = oo and is itself embedded in a region of purely metal- 
lic states. As can be seen from Fig. [HI these metallic 
states occur, for a given fixed interaction strength, only 
for sufficiently large doping. The QMC phase diagram, 
presented here, is in qualitative agreement with the re- 
sults of Riiegg et al. which were obtained based on 
variational methods. As a remark we add that, since, in 
the OSMP at half filling, the wide band is in a non- 
Fermi-liquid state,— by continuity one also expects a 
non-Fermi-liquid state in the OSMP at low doping. Sig- 
natures of such a non-Fcrmi-liquid state away from half 
filling will indeed be found and discussed below. 
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FIG. 7: (Color online) Phase diagram of the two-band model 
for densities 2 < n < 2.6. Within the range of parameter 
values, for which both bands are metallic, the solid green (or 
patterned red) area indicates that the wide (or narrow) band 
contributes stronger to the doping, implying n„ > Un (or 
Un > Tiv,). At the boundary (dashed line), both bands are 
equally filled. 




FIG. 8: (Color online) DOS calculated for various particle 
densities near half filling as a function of frequency (a) for the 
wide band in the two-band model with interaction U = 2.4, 
and (b) for a single-band calculation with identical parameters 
(bandwidth W = 4, interaction U = 2.4, and temperature 
T = 1/40). 



Our phase diagram Fig. [5] is in reasonable general 
agreement with the ground state phase diagram obtained 
by Inaba and Kogai^ in the SU(2) symmetric limit (also 
for J = U/A): the critical interactions at half filling are 
only slightly larger in the latter case (by about 10%); 
however, at strong interactions U ^ 3, the OSM phases 
seem to extend over larger filling ranges for SU(2) sym- 
metric Hund rule couplings. Unfortunately, the coarse 
filling grid employed in Ref. H docs not allow for a de- 
tailed comparison at the scale of Fig. [6l Still, we may 
conclude that the phase boundaries do not sensitively 
depend on the specifics of the Hund rule couplings (with 
or without spin- flip and pair hopping terms). 

The purely metallic phase in Fig. [6] can be subdivided 
into two regimes, depending on which of the two orbitals 
is stronger doped (cf. Fig. [T|), as shown in Fig. [71 Ev- 
idently, the OSMP is entirely embedded in a region in 
which the wide band hosts the larger fraction of the dop- 
ing, which is itself embedded in a parameter region in 
which the narrow band contains the larger doping frac- 
tion. The boundary between both sections of the metallic 
phase is marked by a dashed line, which corresponds to 
the crossing points of the density curves for the narrow 
and the wide band at fixed interaction of Fig. [1] Inter- 
estingly, it is quite easy to predict the large-C/ asymp- 
totics of the (dashed) cross-over line: it should converge 
to n = 3 for U oo. The reason is that equal filling 
of both orbitals is asymptotically expected in all Mott 
phases, in particular also at three quarter filling (with 
Uc ^ 20, cf. Sec. Ill Bl) . Thus, the cross-over line should 
either continue, beyond Uc, as the Mott phase transition 
line or converge towards that line for U > Uc- 



III. DYNAMIC OBSERVABLES AS A 
FUNCTION OF DOPING 

In this section, we discuss dynamic observables: first 
the spectral function (at general frequencies), then, in 
particular, the spectral weight at the Fermi level and 
the functional dependence of the Z-factor, which cor- 
responds to the quasiparticle weight in a Fermi-liquid 
phase but loses this interpretation in the OSMP or an 
insulating state. The spectral function and the spec- 
tral weight at the Fermi level are calculated using the 
maximum entropy method (MEM). Finally, the extent 
of non-Fermi-liquid behavior will be analyzed on the ba- 
sis of imaginary-time self- energies. 



A. Spectral function 

Both in the case of the wide band and for the nar- 
row band a comparison to typical single-band behavior 
is helpful for understanding the doping dependence of the 
spectral function. 

The spectral function of the wide band of the two- 
band Hubbard model, as calculated with the MEM from 
our quantum Monte Carlo data for the single-particle 
Green function, is presented in Fig. [5^ (left panel), where 
it is juxtaposed with the spectral function of a single- 
band model with the same parameters [see Fig. [BJd]; the 
latter corresponds to the wide band in two-band model 
([T|) where the interorbital couplings U' , Jz have been re- 
moved. The spectral function of the wide band in Fig. 
[5^ is characteristic of the typical spectrum of the wide 
band within the OSMP. Minor changes in the interaction 
strength U cause only slight changes in the spectrum. 
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Comparing the result of the wide band of a two-band 
system to the spectral function from a single-band cal- 
culation, one finds that the spectral weight of the wide 
band has a broader distribution. This can be attributed 
to an increased scattering amplitude in the two-band 
model, due to the larger number of interaction channels. 
Moreover, the Kondo resonance is more pronounced in 
the single-band case and the dips between the resonance 
and Hubbard bands are clearly visible in the low-doping 
regime. One also notices a dip in the spectral function 
of the wide band at the Fermi edge in the low-doping 
domain. This phenomenon is well-known from previous 
work2^ and can be attributed to non- Fermi-liquid behav- 
ior. In fact, the same kind of dip is found in the Hubbard- 
III solution or, equivalently, in the exact solution of the 
half-filled Falicov-Kimball model, which describes a non- 
Fermi liquid. Note that the dip persists to about n ~ 2.1, 
i.e., the boundary of the OSM phase (at U = 2.4); be- 
yond that density, it cannot be distinguished from the 
(inevitable) numerical noise in the spectra. We will fur- 
ther discuss the extent of non-FL features in Sec. IIII CI 

We now consider the spectral function of the narrow 
band in a two-band system, which is again juxtaposed 
with the spectrum of a single-band model at comparable 
parameter values. For this purpose, we define scaled in- 
teraction strengths U'-^^ = U/U^^^ and C/^^) = C//C/i^^ for 
the single- and the two-band system, respectively. Here 
Uc^^ and jyi^"* are the critical interactions for the Mott 
transition of the single-band model and the metal-OSMP 
transition of the two-band system, respectively. The nu- 
merical values of Uc^^ and [/^^^ are approximately given 
by: 

U^^'' « 2.35 (single-band model) 

« 2.10 (two-band model) 

The parameter values of the narrow band and the single- 
band system are referred to as "comparable" if C/^^-* « 
f/'^^; then, the systems have about the same "distance" 
from a Mott transition at half filling. In actual QMC cal- 
culations it has to be taken into account that the critical 
interactions Uc^^ and Uc'^^ depend slightly on the Trotter 
discretization. Moreover, in comparing the two bands, 
a slight mismatch of the available U values is unavoid- 
able, since the required Uc values are only approximately 
known and the ratios [/, resulting from the simulations, 
are known only on in general non-identical grids. 

Note that the need of rescaling arises for the inter- 
actions since the relevant parameter is C/ — 1, which is 
strongly affected even by small changes in Uc in the phase 
region of interest. In contrast, due to the small differ- 
ence of only about 10% between J/i^'' and Uc'^\ a rescal- 
ing of frequencies is not necessary for the comparison 
(but would make quantitative comparisons with other 
data more difficult). Similarly, it is not expected that 
temperature variations of the order of 10% would visibly 
alter the spectra; consequently the comparison is made 
at identical (unsealed) temperatures. It should also be 



noted that a rescaling would not have been appropriate 
for Fig. [51 where the wide-band of the two-band system is 
compared with a single-band system of equal bandwidth 
W ^ 4; in that case, only the two-band system is close 
to an MIT (of the narrow band) while the single-band 
system remains metallic at least for U < 4.7. 

Due to the occurrence of a metal-insulator transition, 
the changes in the spectrum of the narrow band, as the 
interaction is varied through U ~ 1, are far more drastic 
than in the case of the wide band. For the narrow band, 
therefore, we discuss results for several interaction ratios 
U , both below and above the metal-insulator transition, 
which are presented in Fig. [S] (panels a, c, e, and g), where 
they are juxtaposed with the "comparable" single-band 
spectra of panels 5, d, /, and h. 

For interaction values U ~ 0.85, clearly below the 
metal-insulator transition(s), we conclude from Fig. [5] 
(panels a and 6) that the spectra of the narrow band and 
the single-band model show only minor differences. The 
Kondo resonance is well-established in both cases, and 
the overall shape of the spectra is very similar. The dips 
in the spectrum between the resonance and the Hubbard 
bands are clearly more pronounced in the single-band 
case, in particular for larger doping. 

For a relative interaction strength U ~ 0.95, still 
slightly below the metal-insulator transition, we note 
from Fig. ^ (panel c and d) that, close to half filling, 
the Kondo resonance of the narrow band in a two-band 
model has considerably smaller weight and width than 
that of the single-band model, much beyond possible scal- 
ing effects. For the narrow band, more spectral weight is 
shifted from the vicinity of the Fermi level to the shoul- 
ders of the Hubbard bands so that the dips between the 
resonance and the Hubbard bands are, again, less pro- 
nounced for the narrow band than for the single-band 
system and the peaks of the Hubbard bands are closer to 
the Fermi level in panel c than in panel d. 

For fJ ~ 1.02 — 1.05, both the narrow band of the two- 
band system and the single band are in a Mott insulating 
state at half filling, as can be seen from the gapped spec- 
tra at the lowest doping levels in panel e and / of Fig. 
[HI As the doping concentration is increased, the spectral 
weight near the Fermi level is larger for the single than for 
the narrow band, in accordance with our previous obser- 
vations concerning pinning of the particle numbers of the 
narrow band at not-too-large doping. In spite of the pin- 
ning of the particle density (rin ~ 1), we note from Fig. 
[5^ that the shape of the Hubbard bands of the narrow- 
band spectrum experiences some considerable changes at 
low doping. At larger doping, we note that the gap be- 
tween resonance and Hubbard bands is much deeper in 
the single-band case than for the narrow band. Moreover, 
the height of the Kondo resonance of the narrow band at 
the largest available doping levels seems to be reduced 
compared to the single-band case due to a larger num- 
ber of available scattering channels. The evolution of the 
overall shapes with increasing doping can also be char- 
acterized in the following way: in the single-band case, 
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FIG. 9: (Color online) DOS calculated for various particle 
densities near half filling as a function of frequency, both for 
the narrow band in the two-band model (left column) and 
for a single-band calculation (right column). Panels within 
one row have approximately matching scaled interactions U. 
In the two-band case (left column), the total filling is varied 
in the range 2 < n < 2.2 in uniform steps of 0.01, with a 
corresponding offset. In the single-band case (right column), 
the filling is varied nonuniformly: here, the offset is chosen so 
that adjacent spectra with the same offset have (roughly) the 
same partial filling (i.e., Un in the two-band model equals n 
of the single-band model). 




2.3 

n (total filling) 

FIG. 10: (Color online) Spectral weight at the Fermi level for 
the two-band model. The thin dashed curves represent the 
spectral weights at the Fermi level according to Luttinger's 
theorem. All other curves are guides to the eyes only. 



the quasiparticle peak grows essentially within the gap; 
in the two-band case, the narrow-band resonance evolves 
from a shoulder of the (upper) Hubbard band. 

For interaction values U ~ 1.11 — 1.14, significantly 
above the metal-insulator transition, one sees from pan- 
els g and h of Fig.[n]that the behavior, already indicated 
in panels e und /, becomes more pronounced. For larger 
doping (n > 2.10) the upper Hubbard band crosses the 
Fermi level, and a Kondo resonance evolves in both the 
narrow and the single band. For the narrow band, the 
resonance is strongly surpressed, however, and the vicin- 
ity of the Fermi level is dominated by the upper Hubbard 
band. Also, there is not a clear gap between the reso- 
nance and the lower Hubbard band of the narrow band 
spectrum. In contrast, the single-band model develops a 
very sharp Kondo resonance which is clearly separated 
from the lower Hubbard band by a (pseudo) gap. 



B. Spectral weight at the Fermi level 

From our MEM-results for the spectral function we 
determine in particular the spectral weight at the Fermi 
level. Our results for this quantity, both for the narrow 
and for the wide band, are presented in Fig. [TOl Within 
DMFT, Miiller-Hartmann^ showed for single-band mod- 
els that, in the ground state, the spectral weight at the 
Fermi level is not renormalizcd by the interaction as a 
consequence of Luttinger's theorem. His proof is also 
applicable to multiband Hamiltonians, diagonal in the 
band- and spin-indices, such as the J^-model ([1]) consid- 
ered here. The reason is that, for such Hamiltonians, the 
multiband Green function, the Weiss field and the self- 
energy are also diagonal in the band- and spin-indices. 
Hence, Luttinger's theorem applies also to the multi- 
band case, if the self-energy has quasiparticle properties: 
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Im E(w) = 0{uj'^) {lli — > 0) on the real axis. Similar argu- 
ments can be found in Ref. |48| . For comparison, we have 
also plotted, in Fig. [TOl the spectral weight at the Fermi 
level as predicted by Luttinger's theorem (dashed lines). 

It is clear from Fig. [10] that the agreement between 
the spectral weights Aw(0) for the wide band and An{0) 
for the narrow band, both calculated at the Fermi level, 
and the values predicted by Luttinger's theorem improves 
with doping. In particular, there are significant devia- 
tions at low doping between the data and the reference 
curves. There are several reasons for this: First of all, 
the narrow band is non-metallic and, hence, not a Fermi 
liquid at low doping (and for U ^ 2.1). Secondly, the 
data were sampled at finite temperature and the refer- 
ence curve is valid only at T = 0. Third, one expects 
that the bands, even when they are metallic, retain some 
non-Fermi-liquid properties at interaction values which, 
at half filling, arc characteristic of an OSMP. Quantita- 
tively, one deduces from Fig.fTOlthat the metallicity (and 
hence the quasiparticle properties) of the wide and the 
narrow band develop in the doping ranges 2.0 ^ n ^ 2.3 
and 2.3 ^ n ^ 2.6, respectively. For larger doping con- 
centrations and interactions U > 2.6, the agreement be- 
tween the data and the reference curves is quite good. 



C. Self-energy and quasiparticle weights 

We briefly discuss the behavior of the quasiparticle 
weight (QPW) Z, as defined by Eq. ([2]), as a function 
of interaction and density. As stressed before, discrete 
estimate ^ oi Z can, strictly speaking, only be in- 
terpreted as a physical "quasiparticle weight" if both 
bands are Fermi liquids. If one now studies the behav- 
ior of the QPW in the two-band Hubbard-model away 
from half filling for interaction values outside the Fermi- 
liquid regime {U ^2), one expects and indeed finds an 
anomaly: The QPW of the narrow band does not de- 
pend monotonically on the filling [Fig. [Tib ] , as it would 
in a Fermi liquid phase, but is increasingly enhanced at 
n w 2.1 with increasing U. To make sure that the behav- 
ior in Fig. [TTb is not a numerical artifact, these results 
were counter-checked with a different QMC solver. 

Having established that the QPW behaves anoma- 
lously outside the Fermi liquid regime, it is interesting to 
search for possible non-Fermi-liquid anomalies in the full 
frequency-dependent self-energies, from which the QPWs 
have been computed according to Eq. The self- 

energies behave smoothly as a function of the (discrete) 
Matsubara frequencies, down to the first Matsubara fre- 
quency, for all metallic solutions. However, already for 
the moderately large interaction value U = 2.4, a clear 
anomaly can be seen in the real part of the self-energy 
[Fig. [Tib]. In the doping region between n = 2.00 and 
n = 2.40 some intermediate solutions show divergent be- 
havior for — > 0, which strongly differs from the smooth 
solutions for densities n — 2.00 and n = 2.40. Moreover, 
T,{iu!n) shows a non-monotonic dependence on the fill- 
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FIG. 11: (Color online) (a) Quasiparticle weight Zn as defined 
in Eq. ([2]| for the narrow band in the two-band model. Zn 
develops a local maximum as a function of filling in the weakly 
doped regime 2.05 ^ w < 2.15 with increasing interaction U. 
There are no further anomalies at stronger doping (see inset). 
(b) Real part of the self-energy calculated for U = 2.4 at 
discrete Matsubara frequencies for the two-band model. The 
self-energy exhibits an anomaly (i.e., a strong enhancement at 
small frequencies) in the intermediate doping range (2.05 ^ 
n < 2.2). The lines are guides to the eyes only. 



ing at the first four Matsubara frequencies wi, uj2, W3, 
and W4. This contrasts strongly with the characteristics 
found in the metallic regime, where Z{iLUn) is monotonic 
as a function of n. Hence the self-energies, too, show 
anomalous behavior, as was found before for the secant 
estimate [see Eq. ^] or Fig. [TTk . Note that the real 
part of the self-energy has to vanish on the imaginary 
axis due to particle-holc symmetry at half filling, even in 
the insulating phase. 

Let us now come back to the question of how closely 
non-Fermi-liquid properties are connected with orbital- 
selective Mott physics (for anisotropic Hund rule cou- 
plings). In Fig. [51 we have already seen that the low- 
frequency dip in the wide-band spectrum, characteris- 
tic of a non-FL, persists at least up to n = 2.1, the 
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FIG. 12: (Color online) (a) Imaginary part of the self-energy 
Im E(icj„) for various total fillings n at [/ = 2.4 and T = 1/40; 
intermediate minima are seen at low frequencies only for n > 
2.15. (b) Results for Im E computed at different temperatures 
for selected densities (at U — 2.4). The curves are hardly 
distinguishable for n — 2.1; for n = 2.2, the low- frequency 
minimum becomes more pronounced with decreasing T. 

boundary of the OSM phase; however, some traces re- 
main up to about 71 ~ 2.15. So it is not clear, at this 
point, whether the OSMP and the non-FL phase are re- 
ally identical. In order to avoid the ambiguities of an- 
alytic continuation, we will now discuss this issue on 
the basis of the imaginary-time self-energy data shown 
in Fig. [T21 Subpanel (a) shows ImI](icL'„) for U = 2.4 
and T = 1/40 (i.e., the same parameters as in Fig. [8]) 
across a range 2.05 < n < 2.25 of densities. Close to 
half filling, at n < 2.1, ImE(ia;) is a monotonic function 
of Lo which clearly tends toward finite non zero values 
(of order 1) for uj ^ 0, i.e., shows characteristic non-FL 
behavior. In contrast, at n = 2.25, far away from half 
filling (and from the OSMP), the self-energy clearly de- 
cays at small frequencies, consistent with the Fermi liquid 
scenario Imi;(iw) oc 0{uj) + 0{T^). In between, the sit- 
uation is less clear: it is impossible to decide, from this 
data alone, whether the apparent finite value to which 
the data extrapolates for to ^ represents thermal or 
intrinsic non-FL effects. Therefore, we have checked the 
impact of variations of the temperature for two selected 
densities. As seen in Fig. [T^] (b), the curves connecting 
the discrete Matsubara data are on top of each other for 
n ~ 2.1. Evidently, the low-temperature physics is de- 
termined by genuine non-FL scattering in this case; tem- 
perature effects are negligible. In contrast, a clear trend 
with decreasing temperature toward a Fermi-liquid form 
is seen for n = 2.2; however, much lower temperatures 
would be needed in order to verify the full restoration of 
FL properties which would be very costly using QMC. 
This is even more true for n = 2.15 for which ImS(itt'„) 
shows only a very shallow minimum. So we have es- 
tablished that Fermi liquid properties are, indeed, re- 
stored for large enough doping; unfortunately our data 



cannot show whether this boundary matches the narrow- 
band OSMT at low temperatures. However, we have 
clearly shown that the wide band is a non-FL throughout 
the orbital-selective Mott phase. Note that the present 
Matsubara-frequency analysis nicely confirms the conclu- 
sions drawn from the wide-band spectra [Fig. 8] and, 
thereby, the reliability of the maximum-entropy proce- 
dure. 



IV. SUMMARY AND CONCLUSION 

We first summarize our results. In this paper we stud- 
ied electronic correlations within two-band model ([T]) 
using a high-precision QMC-DMFT code, both close to 
and far away from half filling. 

We calculated both static and dynamic properties. 
Among the static properties, we focused particularly on 
band-specific fillings and intraorbital double occupancies 
as functions of interaction and doping. The most im- 
portant result for these quantities is their pinning at the 
Mott insulating state for the narrow band. Our results 
for static quantities in a wide range of interaction and 
doping values were summarized in a phase diagram. In 
the section on the dynamical properties we focused on 
the frequency dependence of the spectral function and, in 
particular, on the density of states at the Fermi level. We 
calculated the density of states, both of the wide and of 
the narrow band, and compared the results to those for a 
single-band system. The main finding here is that, in the 
two-band model, the density of states is somewhat broad- 
ened, and that (compared to single-band spectra) weight 
from the Kondo resonance is shifted into the dips be- 
tween the resonance and the Hubbard bands. A physical 
explanation for this is that the interorbital coupling in- 
troduces additional scattering channels for the electrons, 
which reduce quasiparticle lifetimes and, hence, smear 
out the spectrum. 

In this paper, we followed the common practice in the 
study of orbital-selective Mott transitions and focused on 
paramagnetic phases, excluding magnetic states. Still, 
the true thermodynamic equilibrium states of model ([T]) 
may well be magnetic for the parameters of interest. 
In fact, it is known for the half-filled case (in high di- 
mensions) that the Neel temperature at Ud is about six 
times larger than the critical temperature Td ~ 0.02 of 
the OSMT Q so that the OSMTs observed in para- 
magnetic calculations should be hidden inside the an- 
tiferromagnetic equilibrium state. However, studies of 
symmetry-broken states are notoriously difficult away 
from half-filling due to the large variety of possibilities. 
Also, various mechanisms exist which can suppress or- 
dered phases, such as disorder, spin-orbit coupling or 
longer-range hoppingi^ In particular, finite dimension- 
ality can be expected to suppress magnetic states effec- 
tively in planar systems such as Ca2-xSrxRu04. There- 
fore, and since Mott transitions have been observed^i^ 
to persist as narrow crossovers far above critical tempera- 
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tures, results obtained in the paramagnetic phase should 
be qualitatively valid, at least at intermediate tempera- 
tures. 

To conclude, we demonstrated that, for interaction val- 
ues in the OSMP range 2 < [/ ^2.6, the entire doping 
range < n < 4 is subdivided into three different regimes 
by the behavior of the narrow band. At low doping, 
both the particle number and the double occupancy of 
the narrow band are pinned to their Mott insulating val- 
ues, so that, upon increasing the doping, the wide band 
absorbs all additional electrons or holes; in this regime, 
the wide band shows characteristic non-Fermi-liquid be- 
havior. In the intermediate doping regime, i.e., from a 
critical doping concentration ncrit(C^) onward, the nar- 
row band leaves the pinned state, so that both bands 
contribute to the compressibility and become metallic. 
Third, for a nearly filled (or nearly empty) system with 
n > 3.6 (or n < 0.4), the narrow band is effectively 



a band insulator (or effectively empty), so that, also in 
this regime, only the wide band contributes to the com- 
pressibility. In this third regime, correlation effects be- 
come surpressed and quantum fluctuations vanish. Sur- 
prisingly, in the double occupancy this effect can qual- 
itatively already be seen for densities, not too far away 
from half filling (|n — 2| > 0.5). Taken together, our 
findings give a fairly complete and reliable description 
of the physics of two-band model ([T|) with inequivalent 
bands, in which genuine multi-band and orbital-selective 
features are clearly separated from generic correlation ef- 
fects (that also exist in single-band models). 
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